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We prove that second-order (double- loop) chaotic sigma-delta schemes are stable: 
within a certain parameter range, all state variables of the system are guaranteed to 
remain uniformly bounded. To our knowledge this is the first general stability result 
CN for chaotic sigma-delta schemes of order greater than one. Invariably as the amount 

^ of expansion added to the system is increased, the dynamic range of the input must 

q decrease for stability to be guaranteed. We give explicit bounds on this trade-off and 

verify through numerical simulation that these bounds are near-optimal. 
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1 Introduction 



Analog-to-digital conversion is the process of representing a continuous signal by a bitstream, 
or a sequence consisting of +l's and — l's. Sigma-delta quantization schemes often preferred 
in practice for their simplicity and robustness. The standard algorithms however tend to 
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produce periodic output, causing audible idle tones. One approach towards breaking up such 
periodicities is to modify the parameters in the standard sigma-delta scheme and implement 
instead chaotic sigma-delta quantization. However, rigorous results about the stability of 
the double-loop chaotic sigma-delta scheme have remained elusive. It is our goal in this 
paper to prove that under the explicit parameter conditions which we will derive, all state 
variables in chaotic sigma-delta recursions are guaranteed to remain bounded. 



Background on analog-to-digital conversion. Many signals (e.g. sound, light, and 
pressure) are naturally produced as bandlimited signals, being well-approximated locally as 
linear combinations of sines and cosines not exceeding a maximal frequency, or bandwidth. 
As these signals are continuous, they must be sampled and quantized in order to be stored 
in digital format. Recall that a bounded real- valued function / is bandlimited and has 
bandwidth Q if its Fourier transform f(co) = f(t) exp (—2niujt)dt vanishes outside the 
bounded interval [—fl, fi]. The classic sampling theorem gives rigorous bounds on the sam- 
pling rate neessary for accurate reconstruction; namely, it states that a bandlimited function 
can be recovered from its values at evenly-spaced sampling points if the sampling frequency 
is at least twice the bandwidth. For functions / with bandwidth Q = 1/2 (normalized for 
simplicity), the sampling theorem states that for any T > 1, 

where g is any function (a.k.a. low-pass filter) whose Fourier transform is unity inside the 
interval [—1, 1] and zero outside of an interval [— T ,T ] with T < T pQ. 



For a bandlimited function, then, analog-to-digital conversion then amounts to mapping the 
sequence of samples /„ = /(fO to a sequence of discrete values q n in such a way that the 
original function can be reconstructed to within a small error at a later time via 



/(*) = ^ £?«*(*-£) (L2) 



In pulse code modulation schemes, q n is taken to be the binary expansion of f n , truncated 
to its first K bits, sigma-delta (SA) quantization is an alternative quantization procedure, 
popular in practice for its robustness and ease of implementation. For a given input sequence, 
a i^th-order EA scheme dynamically generates a sequence of discrete output in such a way 
that at each step, the quantization error f n — Q n equals the i^th-order difference of a bounded 
state variable v = (v n ) ne z- Informally, this means that the quantization error is pushed to 
high-frequency bands and subsequently suppressed by the low-pass filter g. The following 
proposition provides rigorous bounds on the accuracy of i^th order £A schemes; we refer 
the reader to p] for a detailed proof. 

Proposition 1.1. Suppose f G L 2 (IR) is bounded and has bandwidth Q = 1/2. Let f n = /(^) 



be sampled at rate T > 1, and suppose g in (1.1) is chosen to have absolutely integrable Kth 



derivative \g^ K \t)\dt = H^^Hl 1 < oo. Suppose (f n ) ne z satisfies \v n \ < v max and is such 
that 

f n -q n = A K (v) n (1.3) 
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for some sequence (q n ) n ei, with \q n \ < 1. Here A ( 



E£ (- 1 )'(f)^-^ is the Kth 



order difference operator. Then the quantization error is uniformly bounded by 



fit) 
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The error estimate (1.4) guarantees on the reconstruction accuracy of a Kth order sigma- 
delta scheme. We note that more refined error estimates can be found, for example, in 
[8]. However, in practice the accuracy must be balanced against increased instability and 
implementation costs associated to higher-order schemes. In practice it is common to use as 
low as second-order sigma-delta schemes. 

In the sequel, we will restrict focus to second-order or double-loop SA schemes which gen- 
erate 1-bit quantization sequences q n G { — 1, 1}. For a given input sequence (f n )nez, one 
method for generating output (q n ) n & and state sequence (u n ) n ez to satisfy the second-order 
finite difference equation /„ — q n = A 2 (v ) n is to start with uq — 0, Vq — 0, and iterate for 
n > 1 

q n = Q(u n -i + T^n-l) 

t^n— 1 fn QYi 

V n = Un-x + W„_i + f n ~ q n - (1-5) 

(the symmetric recursion can be implemented for n < —1). Here, 7 is a parameter to be 
specified, and the quantizer Q : R — y {— 1, 1} is the signum function: 



Q(x) 



1, x > 



-1, x < 0, 



Note that when 7 = 1, (1.5) gives the state equations for the typical double-loop SA 
modulator, as in [7] . For the error guarantees in Proposition to hold, we must ensure 
that the state sequence (u n ) ne z remains bounded. This is what is meant by stability of 



a quantization scheme. Treating the system (1.5) as a non-stationary discrete dynamical 
system, Yilmaz proved in [10] that if < (3 < 1, and if 7 is chosen within a certain 
viable /3-dependent range, then stability is ensured and the state sequence can be bounded 
explicitly: \v n \ < Cp. Improved bounds were later provided by Zeng [TT] . 



Chaotic SA quantization. A persistent problem with the SA quantization scheme (1.5), 
and for higher-order schemes as well, is the production of periodic output sequences (q n )nez- 
In audio signal processing, where SA quantization is widely used, such periodicities can 
produce spikes in the frequency spectrum of the reconstructed signal /, and manifest as 
audible idle tones to the listener. Several attempts have been made to modify the standard 



recursion (1.5) to break up periodic output without sacrificing accuracy of the resulting 
reconstruction. One suggested approach is to apply dither, or white noise, to the input f n 
before implementing the recursion [2j. Another approach, not necessarily mutually exclusive 
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to dithering, is to break periodicities by amplifying the state sequence at each iteration, 
considering instead 



u n = A x w n _i + f n -q n 

v n = Xiu n -x + X 2 v n -i + f n - q n ; (1.6) 

This system is called chaotic if either Ai > 1 or A 2 > 1, following several simulation studies 
[13 El IH Ej which indicate the effectiveness of this scheme for breaking up periodicities in the 
output (g n )nez- As we will recall in the Appendix, second-order accuracy |/(t) — f(t)\ < 



C Vmax T~ 2 of the standard scheme is maintained by (1.6) as long as Ai, A 2 are not too large 
compared to the sampling rate T, and as long as the state sequence (v n ) remains bounded. 
While stability for chaotic single-loop SA schemes (K = 1) has been shown in [5], and 



stability for the double-loop chaotic schemes (1.6) was verified in [31 [7] for constant input 



f n = {3 < 1, rigorous stability results for general bounded input |/ n | < /3 < 1 have been 
absent. Stability depends on the choice of parameters Aj, A 2 , 7, and /3; as noted in [3], certain 
parameters lead to instability. 



Contribution of this work. We give an explicit parameter range over which the second- 
order chaotic SA recursion is guaranteed to be stable, |t> n | < v max , independent of the 
sampling rate. We moreover provide explicit bounds on v max - such bounds are crucial 
because quantization inaccuracies arise if state variables exceed certain device-dependent 
limits. 



2 Main results 



Consider the one-parameter family of second-order chaotic 1-bit SA recursions 

q n = Q(u n -i + 7^n-i) 

Un = A« n _i + f n - q n 

V n = Am„_i + V n -! + fn~ <?n, (2.1) 

with quantizer Q returns the sign of its input as before. Then the following stability result 
holds: 

Theorem 1. Fix A > 1 and fix < a < 1 so that 

:= a(l-a)-2y^(l + a)(A-l) > Q 
1 — a + 2y/2(l + a) (A -1) 

Consider a sequence (/ n )nez uniformly bounded by sup n |/ n | < (3 < 1 as input to the chaotic 
double-loop SA recursion (2.1) with multiplier 



If u o = ^o = 0, then the generated state sequence (u n ,v n ) remains bounded and, in particular, 

1-a 



\v n \ < + V for all neZ. 



Note that when A = 1, i.e. we implement the standard recursion, = a and we recover the 



stability results of Yilmaz, (5.1). 




1 1.02 1.04 1.06 1.08 1.1 



Figure 1: Region of stability for the double-loop chaotic recursion (2.1): if the input is 
bounded by (3 and implemented with parameter A > 1, then the state sequence is guaranteed 
to be bounded as long as (A, (3) falls under the curve and 7 is fixed according to Theorem [T] 

While the parameter 7 is assumed fixed in Theorem [TJ we prove a more general stability 
result in Section where 7 can assume any of an interval of values. While the admissible range 
for (3 as a function of A is not derived explicitly in Theorem [T] , we numerically derive this 
range, as plotted in Figure [TJ and see for example that for positivity of (3 in (2.2), necessarily 
A < 1.085. Smaller values of (3 imply tighter bounds on the state sequence. 



3 Numerical experiments 



In this section, we analyze the optimality of the bounds of Theorem [T] via numerical simu- 
lations. 



In Figure |2j we run (2.1) with constant input sequence f n = (3 over a range of values 



of (3, and over a range of values for A > 1. We let 7 be assigned according to Theorem 
[T] The curve (3 b serve d traces the maximal observed value of (3 for which the state 
sequence {v n ) remains bounded. By "bounded", we mean that \v n \ < 1000 for 1 
million iterations, a bound chosen from repeated observations. We compare f3 bserved to 
(^theoretical , the largest value of (3 for which Theorem [l] guarantees stability (in agreement 
with the curve in Figure [T|. Clearly there is a gap between (3 b ser ved and ^theoretical- 
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Figure 2: Theoretical and empirical stability thresholds for the chaotic double- loop SA 
scheme as functions of the expansion A and level of constant input f n = /3. The parameter 
7 = 7(/3, A) is as in Theorem [T] 




1 1.05 1.1 1.15 1.2 1.25 



Figure 3: Empirical stability threshold for the chaotic double-loop SA scheme as a function 
of the expansion A and level ofconstant input f n = /3. The parameter 7 = 1 is fixed. 



The true gap between experiment and theory may be smaller, though, because ^observed 
traces the threshold for stability of constant input only. 

• In Figure [3] we repeat the procedure from Figure |2j but now we fix 7 = 1 and trace the 
minimal observed value of /3 for which the state sequence (v n ) fails to remain bounded. 
We see that, at least for constant sequences f n = (3, the chaotic SA recursion seems 
to be stable for a large range of 7 and over a larger range of A than that guaranteed 
by Theorem [T] 

• The set-up in Figure [4] is analogous to those of Figure [2] and Figure [3j except we now 
compare the theoretical bound on the state sequence v max = sup n \v n \ given by Theorem 
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Figure 4: Theoretical and empirical bounds on the state sequence (v n ) for the chaotic double- 
loop SA scheme as functions of the expansion A and level of constant input f n = (3. 



[TJto the maximal empirical bound on \v n \ at the stability threshold ft 'observed when 7 is 
varied as a function of (3 and A according to Theorem [TJ (open circles). We also plot 
the maximal empirical bound on \v n \ at the stability threshold corresponding to 7 = 1 
(open squares). The theoretical and empirical bounds nearly match, suggesting that 
although our bounds on the parameters (3 and A might be conservative, our theoretical 
bounds on the state variable v max are near-optimal. Note that there is no contradiction 
in observing larger values of v max than the theoretical upper bound given by Theorem 
[TJ as the empirical bound is computed using a larger range for (3. 



4 Set-up and Notation 



To prove Theorem [TJ we first introduce the notation on which the proof relies. We consider 
input sequences (f n )nez that are bounded |/ n | < /3 < 1. We will derive a family of stable 



positively-invariant regions for the chaotic map (2.1 ) as functions of the chaotic factor A and 
a parameter a £ [(3, 1), similar in spirit to the approach of [ID]. We denote by S n := \ f n — q n \ 
the quantization error at time n. Since < /3 < a < 1, this error will be confined to the 
interval [1 — a, 1 + a}. We set 5l = 1 — a and 5h = 1 + «• Consider 

Sf(u,v) = (u-5,u + v-5) (4.1) 

and 

S 5 r (u,v) = (u + S,u + v + S). (4.2) 
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The update rule in (2.1) can be rewritten as 

{Sf n (Xu n - 1 ,v n - 1 ) if q n = 1, 
S^ n (Aw n _i,tf n _i) if q n = -1; 

we will use the shorthand notation 

(u n ,V n ) = S(XUn-i,V n -i). (4.3) 



Let us define the functions 

Bo(u) 



^u 2 + \u + C, u<0, 

2±u 2 + \u + C, u>0. 
^■u 2 + \u-C, u<0, 

^u 2 + ±u-C, u>0, 

where C is a positive constant to be specified. Note that B\ and B 2 are convex and concave, 
respectively; thus the region R = {(u,v) : B 2 {u) < v < Bi(u)} is a convex set, and is 
illustrated in Figure [3j Consider now the action of S from a point (u n ,v n ) in R. Depending 
on whether jv n > —u n or ^v n < —u n , a left move Si or right S r will be applied to determine 
(u n +x,v n +i). This suggests we split R into two subsets, particularly: 

Ri = {(u,v) : v < Bi(u), v>B2(u), > — u] 
R2 = {(u,v) : v < Bi(u), v>B 2 (u), jv < — u] 
R = R^R 2 . (4.4) 



As displayed in Figure [3j the graphs of B 1 and B 2 are denoted by T Bl and T Bi respectively. 
We also illustrate the graph of the line Lq : u + 7t> = 0. The left-most intersection point of 
and Tb 2 is denoted by Po = (-Mo^o), while the intersection of L and Tb 1 is denoted 
by Pi = (— Ui,Vx), and the intersection of Lq and V b 2 by P 2 = (u 2 ,—v 2 ). Because R is 
symmetric about the origin, we have that (u 2 ,v 2 ) = (ui,Vi). A direct calculation yields 
u = {2C5h5l) 1 ^' 2 ■ Another direct calculation gives 

7 = 7 1 2 \ -wr - (4-5) 



5 Proof of stability for the chaotic quantization scheme 



We now prove some supplementary results. It is through these supplementary results that 
we obtain the bounds on parameters (A, 0, 7) that will guarantee the stability of the chaotic 
SA scheme. 
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Figure 5: An invariant set R = R\ U i?2- 



In [10], Yilmaz proved the following stability result for the map (4.3) in the standard setting, 
when A = 1. 



Proposition 5.1 (Yilmaz). Consider the map (4.3) with A = 1. Suppose that \ f n \ < a < 1. 

(5.1) 



Suppose that the constant C in (4.4) satisfies 

2(1 + a 



C > 



1 — a 



and, moreover, that the multiplier 7 lies in the range 



1 + ot ^ ^ J2C{l-a){l + a) - (1 + a) 

< 7 < 



C-(l + a)-'- Ca+ y/C(l -a){l + a)/2 
Then (u f , v') = S(u, v) e R for any (u, v) G R. 



(5.2) 



We observe that C, reflecting the size of the invariant region R, must blow up as a approaches 
1; in practice it is preferable to fix a = a max away from 1 so that state variables remain 
sufficiently small. We will leverage this result to prove stability of the map when A > 1 by 



considering the same family of invariant sets (4.4), but restricting the dynamic range of the 



input to |/ n | < /3 < a, where (5 = /3(A) is diminished to counterbalance the expansion A. 



In particular, we fix a small parameter e > and take (3 



1+s 



; note that [1-0,1 + 0] C 



We first show that under slightly more restrictive conditions than those in Proposition 5.1 
Si(Ri) lies under the graph of B\ even when A > 1. 
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Lemma 5.1. Fix a < 1 and A > 1, and fix e G (0, a). Suppose that C and 7 satisfy (5.1) 
and (5.2), and let 5l = 1 — a and 5h = 1 + ex. Suppose in addition that 



or equivalently, that 
If (u, v) G R and (u', v') = Sf(Xu, v) for 5 G [jztJl, j^6 h ] , then v' < B x {v! 



no < (5.4) 



Proof. Suppose (u,v) G R. Let 5 G [jz^Sl, jt^^h] and consider 5' = 5 — (A — l)u. A 
straightforward computation shows that Sf(Xu,v) = Sf(u,v); that is, a left-iteration by 5 
of the map with expansion A > 1 is equal to a left-iteration by 8' of the standard map with 
A = 1. Since C and 7 are admissible for the standard map by assumption, the stated result 



will follow by Proposition 5.1 as long as 5' G [Sl,Sh]- Observe that because (u,v) G R, 
—Mo < u < uo and v < B\(u). Then 5' = 5 — (A — l)u < 5h for all 5 G [jz^Sl, j^^h] 
and all u G [— uo,uo] if and only if j^Sh + (A — 1)uq < 5h- Rewe warranging, we see that 
this inequality is implied by the stated upper bound on uq. An analogous argument verifies 
further that 5' > 5^, we leave the details to the reader. □ 



Lemma 5.1 gives conditions under which Si(Ri) lies under the graph of B\\ to complete the 
argument that S(R) C R, it remains to show that under a subset of the same conditions, 
Si(R\) lies above the graph of B2. Indeed, by symmetry of the action of S and of the region 
R, Si(R\) C R implies S r (R2) C R. To show that Si(Ri) lies above the graph of B 2 , we 
follow a similar line of argument to that used in [TU] when A = 1. 



Proposition 5.2. Fix cx,C, and e as in Lemma (5.1) and consider the associated invariant 



set R as defined in (4.4). Let 5h = 1 + ct and 5l = 1 — ex. Let (5 = (a — e)/(l + e) and 
suppose that 

Then S t (Ri) C i? /or any 5 G [1 - /3, 1 + 0\ . 



Proof. Lemma (5.1) established that Si(Ri) lies under T^, therefore it is sufficient to show 
that Si(Ri) stays above T B2 . It is easily seen that if v\ > v 2 , then (u,Vi) and (u,v 2 ) get 
mapped to (u',v[) and (u',v' 2 ), with > t>2. Therefore, if we write A = {(u,v) : v = 
B 2 {u) and u 2 <u < Mo}, then we only need to ensure that the map of A, as well as that of 
the line segment connecting P\ to P 2 , stays above Tb 2 - Additionally, for the line segment it 
is sufficient to check just the end points Pi and P 2 because the map Si is affine in (u, v) due 
to the convexity of R. For each end point and for A, we only need to check d~' L := 1 — (3 and 
5' H := 1 + (3 since the map Si is affine with respect to 5. 



10 



1. Case Pi. Because < 5' L < S' H , and because the two points ( — (Ami + S' L ),Vi — (Xui + 
5' L )) and ( — (Ami + S' H ), V\ — (Ami + 5' H )) lie on the same line with positive unit slope, 
whereas B 2 is decreasing for u < 0, both points lie above Tb 2 if the latter point is 



above r^. The condition (5.5) along with the upper bound (5.4) imply: 



Ui < Uq < 



e5r 



< 



eS 



H 



5h - 6' 



H 



(5.6) 



A - 1 - (A - 1)(1 + e) A - 1 
The third inequality follows because e < a. Rearranged, this is Ami + 5' H < U\ + Sh, 



which combined with (5.5) implies that — Uq < —(Xui + 8' H ). 



By construction, Bi(u) — B\{u — Sh) — (u — 5h)- Thus 

vi - Ami - 5' H = Bi{— Mi) - (Ami + S' H ) 
= B 1 (-(u 1 + 5 H )) + (ui + S H ) - (Ami + S' H ) 
> J B 1 (-(n 1 + ^)) >B 2 (-( Ul + S H )) 
>B 2 (-(Xu 1 + 5' H )) 



(5.7) 



the third to last and last inequalities follow respectively from (5.4) and the decrease of 
B 2 when u < 0. 



2. Case Pa: Recall that P 2 = —P± by symmetry of B\ and B 2 . Since B 2 is increasing for 
all m > 0, it is sufficient to show that < Ama — S' H < Ama — S' L < ma for all points Pa- 



The left-most inequality is true because ma > Mi > 5h following (5.5). All we have to 
show is that Ama — S' L < ma, or equivalently, that ma < m < S' L /{X — 1). But this is a 



consequence of (5.4). 



□ 



With Proposition |5.1[ Proposition 5^ and Lemma 5A in hand, we can now prove our main 
result. 

Theorem 2. Fix a < 1 and A > 1 such that A < 1 + 2(1+") ' an< ^ su PP ose that e is in the 
range 

2(l + a)(A-l) 



1 — a 



<e<a. 



Let Sh — 1 + ol and 5l = 1 — ot, and suppose 7 is in the range 



26UX - 1) 



e*S L -25 H (X-iy 



< 7 < 



e_Sj - SjMX 1) 
5 H {X - 1) 



Consider the chaotic double-loop SA scheme (2.1) with bounded input sup n |/ n | < /3 < 1, 
where (3 = If u = v = 0, then the state sequence (u n ,v n ) remains bounded for all 

n GZ; in particular, \v n \ < 2 Sh(\-i) 2 
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Proof. 1. The bound 25h < Uq implied by (5.5) is equivalent to the condition C > 2 f^ L , 



the lower bound (5.1). Combined with the upper bound (5.3), 

25h < r < eHh 

8 L ~ ~ 25 H {\-lf 
The lower bound in this expression is smaller than the upper bound as long as 

2MA-1) 

< e < a, 



(5i 



(5.9) 



and this interval is nonempty only if 



A < 1 + 



a(l — a) 
2(1 + a) 



(5.10) 



2. Given u\ and C, we can express 7 = u\/ B\{— u\) as 

7 



C - (l/2)m - l/(25 H )ul 



The range (5.5) for u\ generates an admissible range for 7: 

(2C5 H 5 L fl 2 - S H 



Sh 
C-5 



H 



< 7 < 



Ca + (C8 H S L /2y/ 2 ' 



(5.1i; 



(5.12) 



The bounds (5.8) and (5.12 



5.1| so we can apply Lemma 



on C and 7 respectively fall within the range of Proposition 



5.1 as well as Proposition 5721 to conclude that if the parameters 



(a, A, e, C, 7) satisfy @, flOp , (|5l0j), and ([511] ), and if < /3 = (a - e)/(l + e), then 
(u',v') = S(u,v) G i? if (w,f) G -R. Theorem |2| follows by maximizing the interval (5.12) for 
7 using the bounds on C, and by noting that max( U] „) efl |u| < C + □ 



Remark 3. Theorem^ is a special case of Theorem 2, obtained by setting e 



2(l+q)(A-l) 
1-a 



6 Discussion 



In this note we proved that second-order chaotic EA schemes (2.1 ) are stable as long as A > 1 
is not too large and the dynamic range ft of the input is sufficiently small. Our stability 
analysis can be extended to a more general setting and is presented in a limited scope for 
the sake of clarity. For instance, our stability analysis also holds for tri-level quantizers in 
(2.1) ; i.e. quantizers Q : IR — > { — 1,0,1}. The motivation for using quantizers having a 
0-output state is to reduce power: in audio processing, over long stretches of input f n ~ 
(such as between songs on a recording), it is desirable that in response q n = so that the 
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analog-to-digital converter can essentially "shut off'. A so-called quiet modification to the 



standard second-order scheme (1.5) was introduced in [9] to have the property that if f n = 



over a sufficiently long stretch of time, then q n = in response. This modification is at the 



same time guaranteed to retain second-order accuracy of the standard scheme (1.4). All of 
the stability analysis of this paper carries over to the quiet setting as well. 



It remains to analyze the full two-parameter family of double- loop chaotic recursions (1.6) 
and to improve the stability theory to better match the empirical bounds (i.e., prove stability 
for a larger range of A). In addition to stability for this model, a rigorous study of the 
dynamics of the quantization output in the region of stability would be of interest. In 
particular, it is still open whether the sequence of quantization output must be truly chaotic 
in the strict mathematical sense. 
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7 Appendix 



In this appendix we show that for Ai and A2 sufficiently small, chaotic recursions of the 



form (1.6) really are second-order, in that they retain the second-order reconstruction accu- 



racy with the sampling rate T (1.4) of the standard double-loop recursion. The following 
proposition is adapted from [ID] . 

and consider 



Proposition 7.1. Let f e L 2 (R), supp(f) G [-1/2,1/2], 
samples (/ n ) n ez = {f( n /T)) neIi with T > 1. Suppose that Ai,A 2 in (2.1) are such that 
Ai, A 2 < 1 + 1/T. Suppose further that \v n \ < v max . Then 



n 
f 



^ Vmax ' Cg • T 



where C g is a constant depending on g, and (q n ) n ez is the output produced from (2.1). 



13 



Proof. Noting that u n = v n — A 2 f n _i, 

/(*)-7?Z)«n5(t- 



< 



T -EA 2 (is('- 

+(Ai - 1)(A 2 - 1) 



T 

(A1A2-I) 



r 



The first term is equivalent to the error of a stable s econ d-order recursion (Ai = A2 = 
1) and is bounded by T'~ 2 ||u|| 0O ||<7( 2 )|| i i by Proposition The second term is similarly 
equivalent to the weighted error of a standard first-order recursion; since [A1A2 — 1| < 2/T it 
is bounded by 2T _2 ||f|| 00 ||p( 1 )||£,i. The final term is a 'zeroth'-order recursion and is bounded 
by T -2 ||u|| oo ||0||£ a . □ 



References 



[1] I. Daubechies and R. DeVore. Reconstructing a bandlimited function from very coarsely 
quantized data: A family of stable sigma-delta modulators of arbitrary order. Ann. of 
Math, 158(2) :679- 710, 2003. 

[2] C. Dunn and M. Sandler. A comparison of dithered and chaotic sigma-delta modulators. 
J. Audio Eng. Soc, 44:227-244, 1996. 

[3] R. Farrell and O. Feely. Bounding the integrator outputs of second-order sigma-delta 
modulators. Circuits and Systems II: Analog and Digital Signal Processing, 45(6) : 691- 
702, 1998. 

[4] O. Feely. Nonlinear dynamics of chaotic double- loop sigma-delta modulation. Circuits 
and systems, 1994. 

[5] O. Feely and L. Chua. Nonlinear dynamics of a class of analog-to-digital converters. 
Int. J. Bifurcation and Chaos, 2:325-240, 1992. 

[6] J.D.Reiss and M.B.Sandler. The benefits of multibit chaotic sigma delta modulation. 
Chaos, 11:377-383, 2001. 

[7] A. Z. Mariam Motamed, Seth Sanders. The double loop sigma delta modulator with 
unstable filter dynamics: stability analysis and tone behavior. IEEE Transactions on 
circuits and systems-II: analog and digital signal processing, 43(8):549-559, 1996. 

[8] T. Nguyen and S. Giinturk. Refined error analysis in second-order sigma=delta modu- 
lation with constant inputs. IEEE Transactions on Information Theory, 50(5), 2004. 



14 



[9] R. Ward. Quiet sigma delta quantization, and global convergence for a class of asym- 
metric piecewise affme maps. Nonlinearity, 23(9), 2010. 

[10] O. Yilmaz. Stability analysis for several sigma-delta methods of coarse quantization of 
bandlimited functions. Constructive Approximation, 18:599-623, 2002. 

[11] S. Zeng. Global Tile Attractor of Second Order Single-Bit EA Modulation. PhD thesis, 
The City University of New York, 2008. 



15 



